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Abstract. The ionization signals in the liquid argon of the ATLAS electromagnetic calorimeter are studied 
in detail using cosmic muons. In particular, the drift time of the ionization electrons is measured and used 
to assess the intrinsic uniformity of the calorimeter gaps and estimate its impact on the constant term 
of the energy resolution. The drift times of electrons in the cells of the second layer of the calorimeter 
are uniform at the level of 1.3 % in the barrel and 2.8 % in the endcaps. This leads to an estimated 
contribution to the constant term of (0.291q'o4) % m the barrel and (O^Iq^) % in the endcaps. The 
same data are used to measure the drift velocity of ionization electrons in liquid argon, which is found to 
be 4.61 ± 0.07 mm/^s at 88.5 K and 1 kV/mm. 
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1 Introduction 

The ATLAS liquid argon (LAr) calorimeter |Tj is com- 
posed of sampling detectors with full azimuthaQ symme- 
try and is housed in one barrel and two endcap cryostats. 
A highly granular electromagnetic (EM) calorimeter with 
accordion-shaped electrodes and lead absorbers covers the 
pseudorapidity range \t)\ < 3.2, and contains a barrel part 
(\r]\ < 1.475) [2! made of two half-barrels joined at r\ = 
and two endcap parts (1.375 < |?7| < 3.2) [3J. Each section 

1 The azimuthal angle <f> is measured in the plane transverse 
to the beam axis. Positive <j> is in the up direction. The pseudo- 
rapidity is defined as rj = — ln(tan(#/2)), where 9 is the polar 
angle from the beam axis. Positive rj is for the proton beam 
circulating anticlockwise. 



is segmented in depth in three layers (denoted as layer 
1,2,3). For |ry| < 1.8, a presampler (PS) [1|3], installed in 
the cryostat in front of the EM calorimeter, provides a 
measurement of the energy lost upstream. 

The EM calorimeter plays a crucial role during the op- 
eration of the LHC, since physics channels involving elec- 
trons and photons in the final state form a crucial part of 
the ATLAS physics program. Achieving the required pre- 
cision and discovery reach places stringent requirements 
on the performance of the calorimeter. The uniformity of 
the calorimeter response over a large acceptance is par- 
ticularly important for the overall resolution. This drives 
several design choices for the calorimeter: lead-liquid ar- 
gon calorimetry provides a good energy resolution and ho- 
mogeneity even in the presence of strong radiation; the 
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Fig. 1. Accordion structure of the barrel. The top figure is 
a view of a small sector of the barrel calorimeter in a plane 
transverse to the LHC beams. Honeycomb spacers, in the liquid 
argon gap, position the electrodes between the lead absorber 
plates. 



accordion geometry (see Figure [T]) avoids readout cracks 
between calorimeter modules, thus also providing good 
uniformity. 

In order to equalize the gains of different calorime- 
ter channels, a calibration procedure involving electronic 
charge injection is used. This is however not sensitive to 
intrinsic characteristics of the ionization gaps in the liquid 
argon system, such as variations in gap sizes and LAr tem- 
perature changes. Such non-uniformities can be measured 
from the ionization signals created by charged particles. 
The calorimeter energy response to this ionization is not 
the best quantity for this purpose, because it requires a 
knowledge of the energy of the incoming particle. How- 
ever the electron drift time in LAr, which can be obtained 
from the signal pulse shape resulting from ionizing parti- 
cles that deposit sufficient energy above the intrinsic noise 
level in a calorimeter cell, is a powerful monitoring tool. 
As explained in Section [2J the drift time is also about four 
times more sensitive to changes in the LAr gap size than 
is the energy response. Cosmic muons have been used to 
this end as part of the calorimeter commissioning before 
the LHC start-up. 

The EM calorimeter installation in the ATLAS cavern 
was completed at the end of 2006. Before LHC start-up, 
the main challenge was to commission the associated elec- 
tronics and automate all of the calibration steps for the full 
173, 312 channels. Cosmic muon data have been taken reg- 
ularly for commissioning purposes since 2006. At the end 
of the summer and during autumn of 2008 stable cosmic 
muon runs were taken with the detector fully operational 
and using various trigger menus. In normal data taking 
only 5 samples around the pulse peak at 25 ns intervals 
are taken, but in order to accurately measure the drift time 



32 samples are needed. The pulse height is also relevant, 
since larger pulses are less affected by electronic noise. A 
summary of the detector performance obtained from cal- 
ibration data, cosmic muons and beam splash events is 
detailed in [5]- 

Measurements of the drift time (T 'drift) in the ATLAS 
EM calorimeter using cosmic muon data are presented in 
this paper. These drift times, which are independent of 
the amplitude of the pulses used for their determination, 
can be compared from one calorimeter region to another, 
and thus allow a measurement of the uniformity of the 
calorimeter. 



2 Ionization signal in the calorimeter 

The current resulting from the passage of a charged parti- 
cle through a liquid argon gap has the typical ionization- 
chamber triangular shape, with a short rise time (smaller 
than 1 ns) which is neglected in the rest of this note, fol- 
lowed by a linear decay for the duration of the maximum 
drift time 



T drift — Wg a p/Vdrift, 



(1) 



Jgap 



is the LAr gap width and Vdrift the electron 



where 

drift velocity [Bj. The ionization current, /, is then mod- 
eled as: 



I(t\ Io,Tdrift) — h I 1 — 



Tdrift 



for < t < T, 



drift 



(2) 



where Iq is the current at t = 0. The peak current ampli- 
tude Iq = p ■ Vdri ft is proportional to the drift velocity and 
to the negative linear charge density p along the direction 
perpendicular to the readout electrode, which varies with 
the lead thickness 0- Since the determination of the en- 
ergy is based on the measurement of Iq, it is crucial to be 
able to precisely evaluate and monitor Vdrift- While the 
LAr gap thickness is mechanically constrained, the drift 
velocity depends on the actual conditions of the detector: 
the LAr temperature and density, and the local high volt- 
age. Uniform response in a calorimeter with constant lead 
thickness requires uniform drift velocity in the gaps. 

At this point it is appropriate to recall that each liquid 
argon electronic cell is built out of several gaps connected 
in parallel: for layers 2 and 3, there are 4(3) double-gaps in 
parallel in the barrel (endcap) respectively; there are four 
times as many gaps per cell in layer 1, given the coarser 
granularity of the readout in the azimuthal direction pQ . 
The parameters measured represent an average of the local 
gaps, both in depth along the cell, and in between the gaps 
forming a cell. 

At the end of the readout chain the triangular signal 
is amplified, shaped and passed through a switched ca- 
pacitor array which samples the signal every 25 ns. The 



2 If the LAr gap increases (as in the endcap) p increases 
slightly on average due to showering in LAr. This is accounted 
for using detector simulation. 
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Fig. 2. Typical single ionization pulse in a cell of layer 2 of the barrel (a) and endcap (b) of the calorimeter. The large red dots 
show the data samples, the small blue dots the prediction and the grey triangles the relative difference (data (S) - prediction 
{g))/ Smax, on the scale shown on the right side of the plot (normalized to the data). 



shaping function (see Section [3]) includes one integration 
and two derivatives. Their net effect is to transform the 
triangular signal in a positive spike, followed by a flat un- 
dershoot, the length of the undershoot being equal to the 
drift time. The net area of the pulse, except for small 
fluctuations due to noise, being equal to 0. Upon Levell 
trigger decision, the samples are then digitized using a 
fast- ADC and recorded [?l[5] . Figure [2] shows two typical 
digitized signal shapes, one for the barrel and the other 
for the endcap. The data samples in each plot correspond 
to a single cosmic muon event in a single cell, and fluctu- 
ations of the amplitude in each sample due to noise can 
be observed. The pulses shown pass the analysis criteria 
described in Section^) The prediction is obtained by mod- 
eling the readout chain as described in Section [3] In the 
barrel section, the nominal gap size is constant (2.09 mm); 
in the endcap the gap size changes with pseudorapidity 
(see Figure [3]), so that at larger values of r] smaller gaps 
lead to a shorter pulse undershoot. 



In the ideal case, an electrode is surrounded by two 
identical gaps, one on each side (see Figure 2|). Any modi- 
fication of one of the gaps by a relative fraction x will break 
the symmetry, leading to two different values of drift time 
Tjji (i = 1,2) (Equations S] and [5]). Figure [3] demonstrates 
this effect by showing the total collected current versus 
time in the case where the electrode is at the nominal po- 
sition (Sg a p = /jm) or shifted by 100 and 200 jim. 
This affects the rise at the end of the pulse (between 450 
and 650 ns on Figure 2(a) for example) which is sensitive 
to changes in the gap size over the charge collection area. 
The variation of the drift time inside the cell arises in part 
from the slight opening of the gaps along the accordion 
folds (see Figure Q}, but the bulk of the effect is caused 
by random or systematic displacements of the electrodes 
away from the gap center. Both effects are parametrized 



by the shift parameter 5„ 



x ■ Wgap- This shift pa- 



rameter is limited to a maximum of 400 fim due to the 
honeycomb filling the gaps, however, some modifications 
of electrical field lines (like edge effects) can contribute to 
local enlargements. 
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Fig. 3. Nominal HV (black dots) and nominal gap width w gap 
(blue triangles) versus rj in the 2nd layer of the EM calorimeter. 




Fig. 4. Schematic view of a LAr gap. The nominal position 
of the readout electrode (dashed line) is exactly equidistant 
from the lead absorbers. Any shift with respect to the nominal 
position (solid line) causes an increase of the gap width on one 
side of the electrode, and a decrease on the other side. 
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Fig. 5. Current as a function of time for a perfect centering 
of the electrode (S gap = (im), a shift of 5 gav = 100 /im and 
Sg a p = 200 fim. 



Beside the gap width, w gap , the model of the signal 
takes into account the electrode shift parameter as well as 
possible variations in high voltage on both sides (neglect- 
ing in a first description the bend parts.). The total signal 
can be expressed as a sum of two triangular signals, one 
for each side of the gap, each described by a drift time Toi 
and peak current • I (i = 1,2). Since the drift veloc- 
ity Vdrift in liquid argon follows, for the range of electric 
fields relevant for this study, a power-law dependence on 
the electric field value [9lll0|. with an exponent denoted 
here by a 



V, 



drift 



= K- 



HV 



Wa 



(3) 



the drift time and peak current fraction are given by: 
T m = T drift (1 - x) 1+a (HV 1 /HV nom )- a , 

fn 



h = 



(1 - x)- a (HVi/HV n , 



T D2 = Tarift (l+x) 1+a (HV 2 /HV n , 

fn 



(l + x)- a (HV 2 /HV n , 



(4) 



(5) 



where T dri f t and f nom (f nom = 1 when the bend parts 
are neglected) are respectively the drift time value and 
the fraction of current corresponding to the nominal high 
voltage HVnom, and HVi corresponds to the actual high 
voltage applied on side i. In the barrel the nominal high 
voltage is 2 kV; in the endcap, the high voltage varies with 
r\ (see Figure [3]) to cope with the varying gap, ensuring in 
principle a constant drift velocity by keeping the electric 
field constant. For the high voltage distribution, electrodes 
are grouped by sectors of Arj x A<f> = 0.2 x 0.2 and for 
redundancy separated supplies are used for each side of 
the electrodes. While in the vast majority of the sectors 
the high voltage has the nominal value, a few of them are 
operated at lower values, to prevent accidental sparking 
or excess noise. 

Both in the barrel and in the endcap, the nominal 
operating field is close to 1 kV/mm. The range of vari- 



ation of x (up to typically 20 %) induces a correspond- 
ing variation of the operating field of ±20 %. In this re- 
duced range, and for a fixed value of the liquid argon tem- 
perature, 88.5 K, the variation of the drift velocity with 
the field is well described [9"IITU] by a power law (Equa- 
tion [3]) . Fitting the data published in [TT] with such law 
gives cti = 0.316 ± 0.030. Additional information was ob- 
tained with our own data comparing a group of sectors 
in the barrel operated at 1600 V, to the bulk operated at 
2000 V. The ratio of the velocity values obtained, taking 
into account small position dependence (see Section [5]), 
gives: a 2 — 0.295 ± 0.020. Considering these two values, 
and given the low sensitivity of our results to the exact 
value of a (see Section [S]) we decided to use a — 0.3 with 
a systematic uncertainty of 11q'q 2 - 

In the accordion geometry, the electric field in the bent 
sections has a lower value than in the flat parts. This leads 
to another contribution to the ionization signal in the form 
of two smaller triangular signals with a longer time con- 
stant and smaller fbend- The sum of the current fractions 
(fnom + fbend) must be equal to 1; the main contribution 
on Figure [5] is related to the drift time in flat sections, 
the tail at large time (t > 500 ns) is due to the larger 
gap width in the bent sections of the accordion. These tri- 
angular shapes are parametrized (neglecting the electrode 
shift effect) by 



Td3 = Tbend {HV\ / HV n 
f, = l^L (HV^HVn 
Td4 = Tbend {HV 2 / HV n 

h = (HV 2 /HV n 



(0) 



(7) 



In the barrel, the Tb en d and fbend contributions per 
layer are estimated using the GEANT4 simulation of a 
uniform charge density in the gap. These values are given 
in Table [T] for layers 1 to 3 (there are no bent sections in 
the presampler). 



Layer 


Tbend ( ns) 


fbend (%) 


Layer 1 


820. 


4.9 


Layer 2 


898. 


7.1 


Layer 3 


941. 


8.5 



Table 1. 

barrel. 



Tbend and fbend values for the different layers in the 



In the endcaps, for practical reasons a different sim- 
ulation was used, MC GAMMA, where 10 GeV electro- 
magnetic showers have been simulated to predict the drift 
time Tdri ft and to estimate Tbend and fbend ■ A photon sim- 
ulation was chosen since the signals relevant to this study 
originate from electromagnetic showers produced by cos- 
mic muons. The simulated photons were generated with a 
flight direction originating from the ATLAS Interaction 
Point. This differs from cosmic muons which cross the 
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(b) MC Tbend for endcap 



Fig. 6. Monte Carlo simulation for (a) Tdrift and (b) Tb £n d versus r\ for the three endcap layers: layer 1 (red triangles), layer 2 
(black dots) and layer 3 (blue squares). 



calorimeter in a quasi- vertical direction. Both Tdrift an d 
Tbend are plotted in Figure [B] as a function of pseudorapid- 
ity for the three layers. These quantities are obtained from 
the distribution of the local drift time where the contribu- 
tions from straight and bent sections of the accordion are 
clearly distinguished. Figure [B] shows that both quantities 
decrease with increasing 77, following the reduction of the 
gap size. The difference observed between the layers is due 
to the depth variation of the gap size: the gap grows con- 
tinuously from layer 1 to 3 due to the projective geometry 
of the cells. The values for layer 2 lie closer to those of 
layer 1. This is explained by the fact that at the energy 
of the simulated showers (10 GeV), the shower maximum 
is closer to layer 1 than to layer 3. The current fraction 
fbend is also estimated from the simulation for every 77 
cell, with values ranging from 5% to 20% depending on 
pseudorapidity. 

3 Prediction of the ionization pulse shape 

The LAr calorimeters are equipped with a calibration sys- 
tem to inject an exponential pulse of precisely known am- 
plitude onto intermediate "mother" boards located on the 
front face (for layer 1) and back face (for layers 2 and 3) 
of the calorimeter. The exponential decay time of these 
calibration signals has been trimmed to mimic the trian- 
gular ionization pulse shape as closely as possible. Since 
the readout path of the calibration signals is identical to 
that of the ionization pulses, the gain and pulse response 
of the electronics can be measured with the calibration 
system over the full range of signal amplitudes and time 
delays. The exponential calibration pulse properties are 
analytically modeled via two parameters T ca u (inverse of 
the exponential slope) and f ste p (relative amplitude of a 
voltage step coming together with the main exponential 
signal) . 

The main ingredient needed for accurate energy and 
time reconstruction in the LAr EM calorimeter is the pre- 
cise knowledge of the ionization signal shape in each read- 
out cell, from which the optimal filtering coefficients jT2] 



are computed. This knowledge of the ionization pulse shape 
is also necessary to accurately equalize the response across 
cells to account for its difference in shape and amplitude 
with respect to the calibration pulse. The difference be- 
tween the two pulses is due to the slightly different shape 
of the induced current (triangle versus exponential) and 
the different injection point for the currents (electrode ver- 
sus mother board). 

The prediction of the ionization pulse shape relies on 
the modeling of each readout cell as a resonant RLC cir- 
cuit (where C corresponds to the cell capacitance, L to the 
inductive path of the ionization signal and R to the con- 
tact resistance between the detector cell and the readout 
line) and on the description of the signal propagation in- 
cluding reflections, amplification and shaping by the read- 
out electronics. 

In the standard ATLAS pulse shape prediction method, 
Response Transformation Method (RTM) [T3] , calibration 
pulses are used to determine the description of the signal 
propagation and the response of the readout electronics, 
as well as the parameters describing the electrical proper- 
ties of the readout cell, (LC and RC) and the calibration 
signal {T caH and / step ). 

A second method has been developed for the EM bar- 
rel, First Principles Method (FPM) [TJ], where the signal 
propagation and the response of the readout electronics 
are analytically described, and the goodness of the ana- 
lytical description is tuned using the measured calibration 
pulses. 

Both methods need, as an input parameter, the value 
of the drift time in each cell, which can be cither inferred 
from the local geometry of the detector along with the 
actual LAr temperature and high voltage, or measured 
from data pulses as described in this work. Details on the 
two methods, which describe the ionization pulse equally 
well, are given below. 
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3.1 RTM Method 

The properties of the signal propagation and of the elec- 
tronic response of the readout of the LAr EM calorimeter 
cells are probed by the calibration system and can be de- 
termined from the measured calibration pulses. The two 
underlying assumptions behind the RTM [13] are that: 

— the ionization pulse (g p hys) can be numerically pre- 
dicted from the corresponding calibration pulse (g ca u) 
by means of time domain convolution with two simple 
functions, parameterizing respectively the shape dif- 
ference between the exponential and triangular cur- 
rents, and their different injection points in the detec- 
tor, see [15] : 

9phys{t) = Qcaliit) 

# c -i f (1 + ST cali )(sT drift - 1 + e - sT <"f') \ 

(. ^^drift (/step + ST C ali) J 

* C ^{l + S 2 LC + sRc} (8) 

where C^ 1 denotes an inverse Laplace transform, with 
s being the variable in the frequency space. The first 
time-domain convolution corrects for the different sig- 
nal shapes through the calibration pulse parameters 
Tcaii and f s tep and the drift time Tdrift, while the 
second convolution accounts for the different injection 
points on the detector cell, modeled as a lumped RLC 
electrical circuit. 

— all parameters (T ca u, fstep, LC, RC) used in the convo- 
lution functions, apart from the drift time, are directly 
extracted from the measured calibration pulses by nu- 
merical analysis [13] . 

3.2 FPM Method 

In the FPM method, the signal generation is based on 
"first principles" of signal propagation [T3] . All the calcu- 
lations are made in the frequency domain, and when the 
signal at the output of the final shaping amplifier is ob- 
tained, it is transformed to the time domain by using a 
fast Fourier transform [15] . 

After generation at the detector level, a signal is prop- 
agated along the signal cable, taking into account its im- 
pedance, propagation velocity, and absorption by the skin 
effect [14] . A small fraction of this signal is reflected at the 
signal cable-feedthrough transition, while the rest is trans- 
mitted. A second reflection takes place at the feedthrough- 
preamplifier transition. In this model, the feedthrough is 
modeled as a single cable section, with its own impedance, 
skin effect absorption constant, propagation velocity and 
length. The preamplifier is described by a complex im- 
pedance, the real part and the imaginary part {Re[ZpA\, 
Iiti[Zpa\) being both functions of the frequency lj. The 
last element of the chain is the CR — RC 2 shaping ampli- 
fier, described by the transfer function: 



F ah (u) = uj ■ r sh /(l + {lu ■ T sh ) 2 f' 2 (9) 

where r s h is the RC time constant of this element. The 
model accounts for both the directly transmitted signal 
and the reflections up to the second order (i.e. two forward- 
backward reflections and two backward-forward reflecti- 
ons). 

Parameters are taken from construction (cable lengths, 
/step and r ca ii, which were measured for all calibration 
boards |16|). from direct measurements channel-by-cha- 
nnel (resonance frequency ujq = \ j\jLC and R) [T7], and 
from measurements on representative samples (Re[ZpA\, 
Im[ZpA], propagation velocity and skin effect constants). 
The signal cable impedance Z$ and the shaper time con- 
stant r s h were left as free parameters and fitted channel- 
by-channel on calibration pulses [2] . The values obtained 
for Zs and T B h came out close to expectations, giving con- 
fidence in the method which describes calibration pulses to 
1% or better. The relative timing of all calibration signals 
was also reproduced with an accuracy of about 1 ns. 

This method was not extended to the EM endcap be- 
cause not all the necessary parameters have been mea- 
sured with the required precision due to a more complex 
geometry. 

4 Description of the data 

Cosmic muon runs from the data-taking period of Sep- 
tember - November 2008 are used in this analysis, cor- 
responding to a period where the LAr data acquisition 
system transmitted and saved 32 samples of the read- 
out signals. The events of interest are those where muons 
lose a substantial fraction of their energy by radiation 
(the energy lost by dE/dx in layer 2 is in average about 
300 MeV L 5:). These events were triggered using calo- 
rimeter trigger towers over the full calorimeter depth, of 
size Aq x A(f> = 0.1 x 0.1 for |r?| < 2.5, 0.2 x 0.2 for 
2.5 < |r?| < 3.2, and up to 0.4 x 0.4 for 3.2 < \r)\ < 4.9. 
Since the data were collected from cosmic muons instead 
of LHC collisions, trigger thresholds were adjusted accord- 
ingly. For technical reasons, only cells which were read- 
out in high gain (LAr readout has three gains with ratio 
~ 100/10/1) are selected for this analysis. This has a very 
small impact on the selected sample as the energy de- 
posits are typically in the high gain range (energies below 
20 GeV). 

Despite the small rate of cosmic muons depositing sig- 
nificant electromagnetic energy, the number of events re- 
corded during the run period ensured sufficient statistics 
for most of the calorimeter regions, with the exception 
of the high-?y region of the endcaps. The pseudorapidity 
range in this study is hence restricted to |7y| < 2.5. 

To minimize distortion of the signal shape, the energy 
deposited in a cell is required to be well above its typi- 
cal noise value. This is particularly important since the 
drift time is obtained on an event-by-event basis. The 
quantity S max is defined as the amplitude of the most 
energetic sample of the data pulse. The minimal required 
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Layer 


S m ax (ADC count) lower limit 


Vnoise (ADC COUIlt) 


F(MeV/ADC count) 




Presampler 


200 


8.0 


7.0 




Layer 1 


500 


8.0 


2.5 


Barrel 


Layer 2 (\rj\ < 0.8) 


160 


5.0 


10.0 




Layer 2 ([r?| > 0.8) 


100 


3.5 


17.0 




Layer 3 


160 


5.0 


7.0 




Layer 1 


500 


7.0 


3.0 


Endcap 


Layer 2 


160 


4.0 


14.0 




Layer 3 


160 


2.0 


7.0 



Table 2. Cut values for the most energetic sample of the data pulse. The approximate electronic cell noise (a n oise) averaged 
over layer and the approximate multiplicative conversion factor from ADC counts to MeV (F) are given as well. 



values for S max are given in Table [5] for the different lay- 
ers of the calorimeter; these values correspond to about 
1 — 2 GeV. The average noise is also quoted, representing 
between 1 and 4% of the minimal value for S max . Unless 
differently stated, all ADC values are pedestal subtracted. 
The difference of thresholds between the \r/\ < 0.8 and 
0.8 < 1 77 1 < 1.4 regions in layer 2 of the barrel is required 
by a difference in gain. To correct for this effect, the nor- 
malized variable SXJ™ is use d for the selection, denned as 
SS = 4.6 • S max for 0.8 < \ V \ < 1.4, and = S max 

everywhere else. An upper limit of 3900 ADC counts for 
Smax P ms pedestal is also required to avoid saturation. 

As a small fraction of the ionization pulses are dis- 
torted and their drift times cannot be determined accu- 
rately, a set of cuts has been defined to select good quality 
pulses: 

— The data should have a negative undershoot in the 
pulse shape. This is ensured by requiring that at least 
5 samples after the peak have a negative amplitude. 

— In order to prevent pulses with too short an under- 
shoot (as can be the case for signals resulting from 
crosstalk for instance), a condition requires that the 
pulse does not contain more than 42 samples around 
ADC counts at the end of the pulse. This condition 
cannot be applied to the endcap where such shapes 
occur due to smaller drift-time values at high pseudo- 
rapidity. 

For a small fraction (6%) of the LAr EM calorimeter 
the high voltage cannot be safely set to the nominal value. 
The cells belonging to these regions are excluded in the 
following. The numbers of pulses per layer after quality 
cuts are given in Table [3] 



5 Extraction of the drift time 

The 32 data samples Si of each calorimeter cell selected 
by the criteria given in Section [4] are fitted using the pulse 
predictions described in Section [31 scaled by an amplitude 
A max and shifted in time by an offset to'- 



9fit(t': A max , to, Tdrijti — A max 
9phys{t', fnom,Tdrift,X, fbendjTtend) for t > to (10) 

Four parameters are left free in this procedure: the 
drift time (Thrift), the associated shift of the electrode 
estimated as 5 gap = x ■ w gap which is in fact only sensitive 
to the absolute value of x when the high voltage is the 
same on both sides of electrodes, the global normalization 
factor A max and the timing adjustment to- The optimal 
set of these four parameters is estimated using the least 
squares method to minimize the quantity: 



Ql = 



1 



{Si — 9fit{ti', A rnax , to, Tdrifti x)) 



y l — l 



(11) 

where n is the total number of data samples used in the 
fit (usually n — 32), N p the number of free parameters 
(N p = 4), and a no i se is given in Table[2j This minimization 
is perform ed us ing the MINUIT package [18j . 

Figure 7(a) presents the variation of Q 2 with S*XJ™ for 
layer 2 of the barrel. An increase of the Q\ value is ob- 
served when S^™ is larger. The same behavior is observed 
in the other calorimeter layers, as expected. In order to be 
able to apply a global selection to the fit quality indepen- 
dently of the data amplitude, a "normalized" Q\, called 
Q 2 , has been used: 





Layer 


# pulses 




Presampler 


20 k 


Barrel 


Layer 1 


43 k 


Layer 2 


331 k 




Layer 3 


79 k 




Layer 1 


13 k 


Endcap 


Layer 2 


45 k 




Layer 3 


18 k 



Table 3. Approximate number of cosmic muon induced pulses 
in each layer after quality cuts. 



Q 2 = 



{Si — 9fit{ti', A rnax , to, Thrift, x)) 



i=l 



+ {kS max )' 



(12) 

where k is chosen such that Q 2 is independent of S max , 
as represented in Figure [7(b) | The denominator in Equa- 
tion [T^] is the quadratic sum of the noise and of the rel- 
ative inaccuracy of the predicted shape. It represents the 
numerator uncertainty. The values of k are given in Ta- 
ble H] for the different layers of barrel (two methods) and 
endcap. 
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Layer 


&fpm in barrel 


feRTM in barrel 


feRTM in endcap 


Presampler 


0.9% 






Layer 1 


1.1% 


0.8% 


0.9% 


Layer 2 


0.8% 


1.0% 


1.4% 


Layer 3 


0.75% 


1.0% 


1.3% 



Table 4. k values for the different methods in the different regions of the EM calorimeter. 
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Fig. 7. (a) Ql versus SXx™ and (b) Q 2 versus Smax in layer 2 
of the barrel. The black points correspond to the mean value. 



For the measurement of the drift time, the last data 
samples corresponding to the end of the pulse are very 
important. It was noticed that for a small fraction of pulses 
(~ 0.6 % for the layer 2) the fit converges successfully but 
the predicted pulse does not succeed in describing the rise 
at the end of the pulse. This implies an incorrect estimate 
of the drift time. To specifically quantify the quality of 
the fit at the end of the pulse, the variable Z\i as t7 has been 
defined, based only on the last 7 samples: 



Jlast7 



32 „ 

E- 

i=26 



(13) 



Large values of |Z\i as t7| single out pulses with erroneous 
fitted drift times. This effect is also observed with a toy 
simulation, and therefore seems to be an intrinsic feature 



of the fitted function, with a large peak followed by a flat 
tail. 

To remove events for which the end of the pulse is 
badly described by the model, a cleaning selection requir- 
ing |Aast7| < 0-15 and Q 2 < 2.5 (3) in the barrel (endcap) 
is imposed. 

An additional set of cuts on the maximum relative 
residual over all samples is applied for presampler cells, 
where pick-up of oscillatory signals was in a few places 
observed (3 % of the pulses): 



\Si — g fit{U)\max / S max < 10%, 

if the residual is small {\Si — gfit{U)\ 
counts), the cut is relaxed to \Si — gfit{U)\ 
20%. 



• mix < 20 ADC 

lax / Smax <C 



After these selections, the fit parameters are examined 
in more detail. Figure [5] presents the distribution of the 
absolute value of the shift parameter, S gap = 
function of the drift time. 



u gap ; 



as a 



The region in Figure [8 (a) with a drift time Thrift com- 
prised between 380 and 550 ns corresponds to the ex- 
pected range for the drift time in the barrel given the 
resolution of the measurement. The low drift time region 
Tdrift < 380 ns of Figure [8(a)] (0.05% of the pulses) is 
dominated by low-amplitude pulses distributed evenly in 
the calorimeter. A closer examination shows that in about 
80% of the cases for the layer 2 barrel, signals in excess of 
Smax = 1500 ADC counts or cells sampled at medium gain 
are found as first neighbors which corroborates a crosstalk 
hypothesis. 

In the region T dri f t > 550 ns of Figure [8(a)] (0.25% 
of the pulses), some pulses are still significantly negative, 
more than 700 ns after the time of signal maximum. A pos- 
sible explanation is that the energy deposit originates from 
a photon emitted along a bent section, thus having an ab- 
normally enhanced fbend contribution. Unfortunately the 
runs taken with 32 samples do not contain information 
from the inner tracker which would have allowed this hy- 
pothesis to be validated by a projectivity study. Aside 
from these extremely large drift time pulses, there is a 
larger class of pulses which are only somewhat longer than 
normal. They are distributed along specific 77 and <fi direc- 
tions: in the transition regions at \rj\ = 0.8 and between 
the two half-barrels at r\ = (see Section 16.1.11) where a 
slight dilution or leakage of the electric field lines yields 
a larger drift time (this is also observed in layer 1 of the 
barrel); in the intermodular regions in <f> in the upper part 
of the detector, where mechanical assembly tolerances al- 
low for a slightly increased gap at the interface between 
modules due to gravity effects (this is not seen in barrel 
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Fig. 8. Absolute value of the shift parameter as a function of the drift time in the barrel (a) and in the endcap (b), for layer 2. 
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Fig. 9. Distribution of the absolute value of the shift parameter in layer 2 of the barrel (a) and endcap (b). 



layers 1 and 3 which are much closer to the mechanical 
fixed points). 

In the endcap, the cloud of points corresponding to 
the expected Thrift is broader than in the barrel, as can 



be seen in Figure 8(b) it ranges from 300 to 600 ns as a 



consequence of the gap size variation with r\ of the endcap 
design. The fact that the dispersion of |<5 ga p| is larger at 
higher values of Thrift is explained as a consequence of 
the larger gap size: the larger the gap width, the larger the 
displacement of the electrode can be. A few events (0.9% of 
the pulses) are observed at very high values of both T^ri ft 
and \S gap \. They are located at low \rj\ where the drift time 



is very large by construction (see Figure 6(a) ). Their pulse 
shape cannot be completely readout using 32 samples, and 
in particular the rise following the undershoot is partially 
absent, which leads to unphysical values of the shift above 
400 ^m. 

A distinctive aspect of the fit, which is clearly visible 
in Figure O is that it yields a peak at \S gap \ = 0. This 
is mainly explained by noise fluctuations. The superpo- 
sition of two triangles of ionization current with unequal 
length due to an electrode shift (see Figure [5]) can only 
lead to a softening of the rise at the end of the pulse, 



compared to the single-triangle case. If, due to noise, the 
rise is steeper than for a single-triangular shape, the fit 
forces Sg a p to 0. In order to improve the statistical signifi- 
cance of high-amplitude signals and minimize the impact 
of noise fluctuations, it has been decided to weight the 
events by (SX!™) 2 - The results in the following sections of 
this note are produced with this weighting factor. 



6 Results in the calorimeter barrel 



Two parallel analyses have been performed for this part 
of the calorimeter using the two pulse shape prediction 
methods described in Section [3] The analyses agree at the 
level of 0.3%, which provides a good check of the robust- 
ness of these results. In this section the measurement of 
the drift time is presented, along with its implications for 
the calorimeter response uniformity and an estimation of 
the electrode shift. 
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6.1 Drift time measurement in pseudorapidity and 
azimuthal angle 

Results in layer 2 are presented first because the statistical 
uncertainties are lower in this layer (see Table [3]) . More 
refined comparisons between the two methods are then 
possible. The following subsection reports on the results 
in the other layers. The presampler is discussed separately 
due to its different structure. 



6.1.1 Layer 2 of the barrel 

Figure [10] presents the drift time Tdri ft extracted from 
the fit as a function of 77. The results of the two methods 
differ by 0.1 ns on average with an RMS of 1.3 ns. The 
full purple line illustrates the prediction from absorber 
thickness measurements made during the calorimeter con- 
struction J . This prediction is based on the fact that the 
mechanical structure of the calorimeter ensures that the 
pitch (with nominal values shown in parentheses) is con- 
stant to within about 5 /im: 



Absorber (2.2 mm) 
+Electrode (0.280 mm) + v 



"gap 



(2.09 mm) 



, (2.09 mm) 



6.66 mm = (2tt/1024) • R, cos t (14) 

where Ri and 6i are the average radius and the local angle 
of the 1024 accordion-shaped absorbers with respect to the 
radial direction. So if the thickness of the absorber varies 
with rj, the gap will also vary in the opposite direction. As 
the drift time Tdrift is directly related to the gap by: 



Tdrift — T£>o(w gap / Wg a p()) 



1+a 



(15) 



a prediction can be derived for the drift time from the 
variations around the nominal gap size (w gap Q = 2.09 mm) 
associated with Too = (Tdrift) = 457.9 ns. 
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Fig. 11. Drift time as a function of <f> in layer 2 of the barrel: 
using the RTM method (open dots), the FPM method (red 
triangles) and the prediction described in the text (purple line). 



The agreement between the prediction coming from 
precision mechanical probe measurements of the absorber 
thickness and the data is rather good, except in the tran- 
sition regions around 77 = 0, ±0.8 and —1.4, where the 
lower held induces a larger Tdrift- To quantify the agree- 
ment between the drift time measurements from the fit 
and the estimate from the measurement of the absorbers, 
the RMS of the difference between the data points and 
the prediction is computed. This yields a value of 2.9 ns, 
as compared to an RMS deviation with respect to a con- 
stant value of 3.7 ns, excluding the data points around 
the transition region in each case. Comparing bin by bin 
the drift times obtained (Figure ITOl) for the negative and 
positive values of 77, one gets a distribution with a mean 
of 3.4 ± 0.2 ns and RMS of 1.7 ns. The predicted value is 
1.5 ± 0.2 ns. 

The Tdri ft distribution as a function of <j) is presented 
in Figure [fl] for both methods. There is a small differ- 
ence between the < ((456.8 ± 0.3) ns) and > 
((458.3 ± 0.3) ns) regions: a (0.3 ± 0.1)% relative effect 
consistent with sagging and pear shape deformation of 
the calorimeter. No significant variations are observed in 
the absorber thickness measurements. The distribution of 
the results is also rather uniform when looking at the two 
half-barrels separately. The RMS of the </> distribution is 
smaller (1.8 ns) when the two half-barrels are combined, 
than for the r; < (2.8 ns) and ij > (3.1 ns) half-barrels 
separately. This may be due to the existence of small 4> 
modulations with opposite phases in the two half-barrels 
that appear to be more visible in layer 3 (see Figure [T^]) . 



6.1.2 Other layers of the barrel 



Fig. 10. Drift time as a function of r\ in layer 2 of the barrel: 
using the RTM method (open dots), the FPM method (red 
triangles) and the prediction described in the text (purple line). 



The distribution of Tdrift as a function of 77 for layer 1 
is displayed in Figure [T3] The results of the two methods 
differ by 1.3 ns on average, with an RMS of 4 ns, and at 
some points by up to 7 ns. The front layer is particularly 
intricate because of the large relative variations of the cell 
depths which present a discontinuity at |?7| = 0.8, inducing 
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a corresponding variation of the cell capacitance and bent- 
to-flat ratio. Given that the two methods differ in their 
estimation of the cell capacitance, such a difference is not 
unexpected. 

In Figure Q21 a drift time modulation with cf> is clearly 
visible for \rj\ < 0.5 (and equally present in both meth- 
ods) in both half barrels of layer 3. While the source of 
the modulation is so far unexplained, the fact that the 
modulations in the two half-barrels are opposite in phase 
is expected, since one of the half-barrels was rotated by 
180 degrees about the vertical direction for final integra- 
tion. 



of the calorimeter. In addition, this gap varies with rj; the 
values for the 4 regions are given in Table [5j The effect on 
the fitted drift time can be immediately seen in Figure 
The prediction superimposed on the measured distribu- 
tion is normalized to the region 0.8 < \rj\ < 1.2. Good 
agreement within 1% is observed between the measured 
and expected drift times as a function of 77. As there are 
no bent sections in the presampler, the pulse description 
is simpler than in the case of the other layers. While the 
variations in r] are large, the <p dependence of the drift 
time is negligible. 
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Fig. 12. 2D map of Tdrift in (?7,0) f° r layer 3. The empty bins 
correspond to sectors with non nominal HV. 
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Table 5. Gap values in presampler. 
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Fig. 14. Drift time as a function of 77 in the presampler barrel 
using the FPM method (red triangles). The full purple line 
represents the prediction normalized to the region 0.8 < \r)\ < 
1.2, using Equation 1151 and the gap values given in Table [S] 
The empty bins correspond to sectors with non nominal HV. 
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Fig. 13. Drift time as a function of r\ in layer 1 of the barrel: 
using the RTM method (open dots) and the FPM method (red 
triangles) . 



6.1.3 Presampler 

The presampler is constructed differently from the other 
layers of the calorimeter. It is made of narrow flat elec- 
trodes. The size of the gaps is slightly smaller than else- 
where, leading to values of Tdrift lower than in the rest 



6.2 Response uniformity 

The reconstructed value of the energy deposited in the cal- 
orimeter by an electron or photon should be independent 
of the position of its impact on the calorimeter. The non- 
uniformity coming from local variations of the response 
due to gap fluctuations can be determined using the drift 
time measurements. This study is done only for layer 2, 
which is the main contributor to the energy response of the 
detector as it collects ~ 70% of the total electromagnetic 
signal in the calorimeter. 

Figure fTS] shows the distribution of the drift time av- 
eraged over groups of 4 x 4 cells corresponding to an area 



2.3 



of 0.1 x 0.1 in Ar\ x Ac/> plane. This area represent a typ- 
ical transverse size of a single particle shower. The aver- 
age of the statistical uncertainties on Tdrift obtained for 
pulses within the various 4x4 groups is 1.25 ns, well be- 
low the dispersion of the determined Tdri ft values of the 
groups (the RMS is 5.85 ns). From the measurement of 
drift times, the systematic dispersion of the gaps can be 
estimated and its impact on the calorimeter energy re- 
sponse can be assessed. 
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Fig. 15. Drift time uniformity between groups of 4 x 4 cells 
(Arj x A<f> = 0.1 x 0.1) for barrel layer 2. 



The drift time uniformity, corresponding to the ratio of 
the RMS and the mean value of the local Tdri ft distribu- 
tion amounts to (5.85±0.14)/457.8 = (1.28±0.03)%. From 
the relation between the drift time and the drift velocity 
(Equation [lj, the latter being proportional to the energy 
response, together with Equation [T5l it follows, that the 
drift time uniformity leads to a dispersion of the response 
due to the gap variations of (1.28 ±0.03) %• (a/(l + a)) = 
(0.29 ± 0.01)%. Excluding transition regions in r\ and in 
<fi, the gap variations amount to 5.7/457.4 = 1.25% and 
the impact on the response is 0.28%. Taking into account 
small variations observed in the result when changing the 
weighting, the fit strategy (see Section [9]) or the pulse 
reconstruction method, a systematic error of 0.03% is ob- 
tained. The uncertainty on a (see Section [5]), treated here 
as external parameter, contributes with a systematic un- 
certainty of 02- Grouping all errors together in quadra- 
ture gives as the final result: (0.29±g g§) %. 



6.3 Electrode shift 

As presented in Section [3J there is some freedom for the 
electrodes to be displaced with respect to their nominal 
positions equidistant between two neighboring absorbers. 
This displacement is expected to be less than 400 /jm ex- 
cept perhaps in the transition regions between modules. 

The electrode shift is left as a free parameter in the fit 
to the data, which yields one value per calorimeter cell. 
Only the average of the absolute value of the displacement 
can be observed. Since a cell consists of several electrodes, 



an effective value is obtained which is a combination of 
the individual movements of each electrode within a cell. 

The local average value for the shift parameter per 
bin of 0.1 x 0.1 is shown in Figure [TBI for layer 2. It indi- 
cates that the bottom half (negative </>) of the negative- 
•q half-barrel has shift parameter values somewhat lower 
than average. Similarly the module azimuthally located 
between 47r/16 and 57r/16 in the r\ > half-barrel presents 
lower shift values. These variations given their distribution 
throughout the detector, are likely to be due to mechanical 
construction issues. 

The shift parameter also covers local variations of the 
"double-gap" within a cell, for example, by the slight open- 
ing of gaps along an accordion fold. This latter variation 
is in general much smaller than the off-centering of elec- 
trodes between absorbers. 
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Fig. 16. (r),4>) map in which \5 gap \ is plotted per bin of 0.1x0.1. 



Smaller values of the shift parameter are expected for 
the presampler compared to the accordion layers, due to 
mechanical constraints on the electrodes which are indi- 
vidually glued in between two precision structural frames 
[2] . The mean value of the shift in the presampler is found 
to be (|<5 9 op|) = 66.5 /im, as compared to 146 /im in the 
accordion section. 



7 Results in the calorimeter endcap 

As was done for the barrel, the endcap results are grouped 
in three different parts: drift time measurements, calorim- 
eter response uniformity and electrode shift determina- 
tion. 



7.1 Drift time measurement in pseudorapidity and 
azimuthal angle 

The drift time Tdri ft averaged over <p is studied as a func- 
tion of -q for each of the three layers of the endcap (see 
Figure [T7|). The two endcaps, A (rj > 0) and C (rj < 0), 
are combined in the figure. A general decrease of Tdrift 
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with increasing pseudorapidity is observed, as expected 
from the corresponding reduction of the design gap size. 
Fewer fluctuations are observed in layer 2, which offers a 
larger cross section to cosmic muon-induced electromag- 
netic showers. In all layers regular steps are observed, 
corresponding to the locations of the boundaries between 
high voltage regions. 
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Fig. 17. Drift time versus pseudorapidity for layer 1 (a), layer 2 
(b), and layer 3 (c) cells of the endcap. Black points are the data 
and red triangles Monte Carlo predictions for photons. The 
vertical dashed lines show the boundaries between different 
high voltage regions. 



The data are compared to the Monte Carlo calculation 
described in Section[2j Good agreement is observed at high 
r], however the Monte Carlo is slightly above the data at 
low values of T) (~ 1— 3 %), which is a more difficult region 
to simulate. 

In Figure 1181 for a comparison, the data points from 
the three distributions of Figure [17] are super- imposed on 
the same plot. An increase of the drift time with the cell 
gap size at fixed rj is clearly observed, with Tdri ft being 
smallest for laye r 1 and highest for layer 3 (see Section [2] 
and Figure 6(a) ). The drift time for layer 2 lies half way 



between layers 1 and 3 in contrast to the Monte Carlo sim- 
ulation (Figure 6(a) ) where the values for layer 2 are closer 
to the values of the layer 1. This difference reflects the 
fact that cosmic muons are randomly distributed within 
the depth of layer 2, while the photons of the simulation 
develop there shower closer to layer 1. 
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Fig. 18. Drift time versus pseudorapidity for the three layers 
of the endcap: layer 1 (red triangles), layer 2 (black dots), layer 
3 (blue squares). The vertical dashed lines show the boundaries 
between different high voltage regions. 

Figure [19] shows the drift time Tdri ft as a function of 
azimuthal angle for layer 2 for the two endcaps. The val- 
ues of Tdri ft for each given pseudorapidity bin have been 
normalized to the average in order to mask the depen- 
dence on 77. Vertical dashed lines indicate the boundaries 
between modules. An asymmetry is visible on Figure 1191 
between positive and negative values of <p: Tdrifti^ > 0) is 
larger (0.996 ± 0.002) than T drift ((j> < 0) (0.980 ± 0.002). 
Since cj) < is the lower half of the calorimeter, we asso- 
ciate this effect to the greater gravitational compression 
of this part leading to slightly smaller gaps than in the 
upper half <fi > 0. 



7.2 Response uniformity 

An estimate of the intrinsic uniformity of the endcap can 
be made in a similar manner as presented for the barrel 
in Section 16.21 The average drift time across a region of 
size 0.1 x 0.1 on the (r],<f)) plane is computed, with special 
care to take into account the varying gap thickness. 

Figure [2U1 represents the distribution of Tdrift/ <T > 
for layer 2. The normalization < Tq > corresponds to the 
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Fig. 19. Drift time normalized to the average value versus cj> for layer 2 of the 7] > (a) and r\ < (b) endcap wheels. The 
black dots are the average per <f> bin and the vertical dashed lines show the boundaries between different modules. 



value (per r\ cell) predicted from a first order polynomial 
fit to the data Thrift in each high voltage region. This 
normalization cancels out the change of the drift time due 
to the nominal design gap size variation with r\. The study 
is carried out only for layer 2 since it contains most of the 
shower energy of typical LHC electrons and photons. In 
addition, more events have been recorded in layer 2 than 
in the other layers, which increases the statistical accuracy 
of the measurement. 

The drift time uniformity of the Thrift (0.1 x 0.1) dis- 
tribution has an RMS of (2.8 ± 0.1)%. To get the pure 
systematic non-uniformity between the 0.1 x 0.1 cells, the 
dispersion within the 0.1 x 0.1 cells, which in this case is 
not negligible, (1.5 ± 0.1) %, is quadratically subtracted. 
These numbers translate to a uniformity of the endcap cal- 
orimeter response due to intrinsic gap variations of (0.54± 
0.02)%. Systematic effects as discussed in Section ® and 
the uncertainty on a (see Section^ increase the error to 
(0.54t°;jS) %. 
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Fig. 20. Drift time uniformity between groups of 4 x 4 cells 
(Arj x A(j> — 0.1 x 0.1) for endcap layer 2. The normalization 
< To > is obtained as a fit to the data using a first order 
polynomial in each HV region to cancel out the influence of 
the gap variation with rj. 



7.3 Electrode shift 

The distribution of the electrode shift as a function of the 
azimuthal angle is presented in Figure [21] for layer 2. A 
rather flat behavior is observed. Vertical dashed lines cor- 
respond to the boundaries between consecutive modules. 
With a finer binning no particular increase of the shift 
is observed at these transitions, even when extending the 
scale to 1000 /im. The average of about 146m is indepen- 
dent of the layer. 
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Fig. 21. Electrode shift as function of cj> for layer 2 of the 
endcap. The black dots are the average per <f> bin and the 
vertical dashed lines show the boundaries between different 
modules. 



8 Drift time and velocity measurements 

To quantify the consistency of the drift time measure- 
ments, the drift velocity (V^ri/t) is studied more closely. 
The drift velocity can be extracted from drift time mea- 
surements if the local gap values are accurately known 
(see Equation [T]) . Both w gap and Tdri ft are designed to 
be constant for the barrel, but varying with pseudorapid- 
ity for the endcap. The variation of the drift time Tdri ft 
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(see Figure 22(a)| does not compensate for the variation 
of Wgap because Thrift ~ w gap ■ ^ n addition, the different 
high voltage regions in the endcap introduce steps in the 
behavior of the drift velocity as a function of rj. 

In order to compare accurately the drift velocities be- 
tween barrel and endcap and for each calorimeter layer, 
they are scaled to a reference field of 1 kV/mm: 



V drift (lkV/i 



2000 V- w a 



T drift \ HVnom ' 2 mm 



(16) 



where HV nom is the nominal high voltage value, w gap is 
taken from the design value and a is the exponent intro- 
duced in Section [2] Figure 22(b) shows the drift velocity 
at the same field 1 kV/mm for layer 2 of the entire calo- 
rimeter as a function of rj. As expected, a rather constant 
behavior is observed over the entire calorimeter. The devi- 
ations from a perfect horizontal line is explained by local 
non-uniformities. Deviations are observed at the transi- 
tion regions at 77=0 and \rj\ = 0.8 and in the crack region 
between barrel and endcap at \rj\ — 1.4, where the field is 
lower. 



rel (88.5 K) and endcap C (88.4 K). This can explain the 
larger drift velocity measured in endcap C (7 7 < 0) with 
respect to endcap A, by ~ 0.6% (see Figure 22(b)), the 
expected difference being approximately 0.5%. 

Figure [531 shows the comparison of Vd r » ft for the dif- 
ferent layers of the barrel and endcaps. The mean values 
of the distributions are also quoted. The errors on these 
means, given the large number of pulses averaged and the 
random nature of the noise dominating the error on single 
measurements, are much smaller than the systematic un- 
certainties (see Section According to Equation [TBI the 
uncertainty in the drift velocity depends on uncertainties 
in both the gap size and the drift time. The former can be 
extracted from an azimuthal and pseudorapidity unifor- 
mity study, giving values smaller or equal to 1% and 2% 
for the barrel and endcap respectively. The latter receives 
contributions from several sources (see Section [§]). The 
mean values of the drift velocity for the different layers of 
the barrel and endcap are given in Table [SI They are all 
compatible within errors, although the barrel presamplcr 
is somewhat below the average. 
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Fig. 22. (a) Drift time and (b) Drift velocity (at E = 
1 kV/mm) versus rj in layer 2. The black dots are the aver- 
age per rj bin. 

The temperature in the endcap A (rj > 0) is slightly 
higher (by about 0.3 K) than the temperatures of the bar- 
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Fig. 23. Drift velocity distribution for the barrel (a) and end- 
cap (b). 



These results can be compared with the measurements 
from [TTJ which give (4.65 ± 0.12) mm/ /is for a LAr tem- 
perature of 88.5 K and provides good agreement with the 
present measurement. 



27 





Layer 


Drift velocity (in mm/ /j,s at 1 kV/mm) 


Barrel 


Presampler 
Layer 1 
Layer 2 
Layer 3 


4.52 ± 0.001 (stat) 

4.62 ± 0.003 (stat) 

4.63 ± 0.002 (stat) 
4.59 ± 0.002 (stat) 


to.or (^t) 

+ -° °A (syst) 

±2:3 (syst) 
+ -°oit (syst) 


Endcap 


Layer 1 
Layer 2 
Layer 3 


4.65 ± 0.002 (stat) 
4.69 ± 0.001 (stat) 
4.59 ± 0.002 (stat) 


(syst) 

±° :i4 (syst) 

i° :3 (syst) 



Table 6. Drift velocity at E — 1 kV/ mm in the different layers of the calorimeter. 



9 Systematic uncertainties 

The different sources of systematic uncertainties affecting 
the measurement of the drift time which have been studied 
are discussed below. The resulting systematic uncertain- 
ties on the velocity are given in Tabled and in Sections l6.2l 
and 17.21 for what concerns the uniformity of response. 



9.1 Comparison of the results obtained in the barrel 
with the two prediction methods 

Two pulse shape prediction methods have been used for 
the barrel. Their results are compared to give an estimate 
of the systematic uncertainty on the prediction. For lay- 
ers 2 and 3, the mean value of the difference between the 
predicted distributions is ~ 0.2 ns and the RMS in the rj 
direction is ~ 1.2 ns which is of the order of the precision 
of the measurement for both methods: hence no signifi- 
cant difference is observed for these layers. For layer 1, 
which also suffers from low statistics, the mean value of 
the difference (1.3 ns) (see Section l6.1.2l) is taken as an es- 
timate of the systematic uncertainty associated with the 
prediction. 



9.2 Different fit strategies 

In addition to the fit procedure described in this paper, 
another approach was also followed in layer 2 of the barrel: 
the cell-based fit. The method consists of fitting simulta- 
neously all the (N) pulses collected in a given cell, using 
a single value for each of the drift time and the shift pa- 
rameter, and N global normalization factors and timing 
adjustments (one of each per pulse). This yields results 
that are similar but not identical to those obtained from 
a weighted average of the individual fits with the weight 
(Smax) 2, F° r instance the average drift time in the case of 
the cell-based fit is 1.2 ns (i.e. 0.3%) lower due to a some- 
what reduced effect of pulses with large Thrift- With the 
cell-based method, which has more statistical power for a 
given fit, the spike at zero visible in Figure [H] is very much 
reduced, confirming that it originates from statistical fluc- 
tuations of the noise leading to a rising slope around 550 ns 
steeper than for a single triangle. 



9.3 Variation of parameters of the cell 

The effect of the uncertainty on the capacitance in layer 2 
of the barrel on the FPM determination of the drift time 
is studied as follows: the capacitance is varied by ±5% 
based on measurements, and a new set of the parameters 
T s h and Zs (defined in Section [3]) are recalculated from 
the FPM calibration fits and used in new fits of the cos- 
mic muon data. A small change in the overall drift time 
scale is observed, but no significant variation in the drift 
time dependence on rj. It should nevertheless be noted 
that when varying the capacitance in either direction, the 
drift times increased by about 3 ns. As discussed in |14j . 
an increase (decrease) in the value of the capacitance is 
partially compensated by a smaller (larger) value of the 
shaper time constant T s h, which leads to only minor vari- 
ations in the pulse shape. 

For the RTM method, the estimated uncertainty for 
the determination of LC and r ca n is less than 3%. The 
T ca i; uncertainty induces an uncertainty of about 0.5% on 
Tdrift, with an additional contribution of less than 0.1% 
coming from the LC uncertainty. 

9.4 Effect of electron attachment 

In the presence of impurities in the LAr medium, drift 
electrons may attach to the impurities with an associated 
lifetime Tu ve , and the signal shape is no longer triangular 
but has the form: 

I(t) = -^e-^'- (l - for < t < T drift 

J- drift \ 1 drift J 

(17) 

Using the Fourier transform of lit), the pulse shape 
is derived by convolution of the various factors affecting 
pulse formation and propagation (see |14j for the general 
case). The data are then fitted with the additional param- 
eter Tn ve . 

Although this new parametrization allows to reduce 
the size of residuals, the values obtained for Tu ve have a 
large dispersion (about 6 jus for both the average and RMS 
of the distribution). Another weak point of this description 
is that the effect is totally absent in the presampler, which 
is in the same liquid bath as the calorimeter. 

A systematic uncertainty of ^J' 5 % in the drift time is 
conservatively estimated from the difference between the 
cases of including or not the Tu ve parameter. The rj de- 
pendence of Tdrift remains essentially the same in both 
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cases. The drift velocity remains unchanged in the pre- 
sampler, but is reduced by 1.5 to 2% in the other layers, 
which would make the presampler and the rest of the bar- 
rel more compatible. 

While this study was made only in the barrel, the esti- 
mated systematic uncertainty is also used for the endcap. 



9.5 Variation of the bent triangle contribution 

The amount of energy deposited in the bent sections of the 
calorimeter is estimated using the simulation. To account 
for possible differences between data and simulation, a 
systematic uncertainty related to the estimate of the frac- 
tion of signal collected in the bent sections is assessed by 
varying the contribution of the triangle associated with 
the bends fbend by ± 20 % based on Table [TJ This test 
was done in a limited section of layer 2 of the barrel. The 
resulting systematic variation of the drift time is =f3 ns, 
as if Tdri ft were compensating the absence of the bent tri- 
angles. It should be noted that the variations of the drift 
time with the relative amplitude of the third and fourth 
triangle (see Equations [4] to [7]) are constant throughout 
the detector; uncertainties on the contribution from bent 
sections should therefore not affect the estimate of the 
intrinsic uniformity, except in layer 1 (see Section [6 . 1 . 2|) . 

The procedure to estimate T Den d and fbend m the end- 
cap requires that the contributions from the bent and 
straight parts of the accordion can be separated using 
the local drift time distribution of simulated 10 GeV pho- 
ton showers. The uncertainty induced by this procedure 
is propagated to the final Thrift value, leading to a 0.2% 
variation that is compatible with the precision of the mea- 
surement. 



9.6 Variation of the parameter a. 

The effect of the uncertainty on the exponent a in the de- 
termination of the drift time in the endcap was studied by 
varying a in the range from 0.30 to 0.39 larger than the 
range determined in Section [2] (0.28 to 0.34). This larger 
range was initially motivated by a previous measurement 
of this exponent during the beam test of the endcap pro- 
totype using 120 GeV electrons, where a value of 0.39 
seemed to describe the data better, however over a larger 
electrical field range than relevant here. The effect of this 
difference (0.30 to 0.39) on the drift time is approximately 
1 ns or about 0.2%, again at the level of precision of the 
measurement. The effect on the retained range (0.28 to 
0.34) would be even smaller. The arithmetic effect of the 
uncertainty on a on the uniformity was considered in Sec- 
tions [O and O 



9.7 Summary of the systematic uncertainties 

The systematic uncertainties discussed above apply to the 
drift time measurements and can be translated into drift 



velocity through the Equation [TJ The drift velocities for 
each layer are summarized in Table [6] 

Averaging over the presampler and layer 2 (barrel and 
endcap) values, for which most of the systematics are un- 
corrected, gives as the final result for the reference field 
of 1 kV/mm and a temperature T = 88.5 K: 



V, 



ref 



(4.61 ±0.07) mm/yits 



(18) 



10 Direct determination of local gap and drift 
velocity at operating point 

Taking advantage of the studies presented above, a some- 
what more global treatment of the data is presented below, 
which allows: 

— to unify the comparison of the local measured gaps, 
and their reference value from construction in both 
the barrel and the endcaps. 

— to obtain for the whole calorimeter the values of the 
drift velocity at the local operating points. 

If the drift velocity were to be fully saturated, i.e. in- 
dependent of the electric field, a measurement of the drift 
time would trivially give the associated local gap using 
Equation [TJ In the situation analyzed here, the drift ve- 
locity depends weakly on the electric field, with a power 
law already given in Section [3] (see Equation [3J . 

Using Equations [TJ and [3] rewritten below as 



Vdrift — V re f 



HV WggpQ 



»' 



(19) 



it is possible to express both the local velocity and the 
local gap, as functions of the measured Thrift- 



Wo 



[A ■ Tdrift 



l/(l+a) 



Vdrift — 



t: 



a/(l+a) 



(20) 



(21) 



drift 



with A = V re f 

below uses 



HV 



HV \ ' W gapO 

a = 0.3, w 9 ap0 = 2 mm 
normalizes the drift velocity at 1 kV/mm to the average 
V re f = 4.61 mm/ /is, as reported in Section [5] The effect 
of the shift (x ~ 0.1) was estimated to bias the above 
analysis by less than 0.2 % on the extracted gap value, 
and is therefore neglected. Data for the endcaps have been 
corrected for the temperature difference, and rescaled to 
88.5 K. 

The additional information yielded by this analysis 
shows directly how the ratio of the measured gap to the 
designed gap varies as a function of position in the detec- 
tor. Figure [M] shows the relative difference between the 
calculated and design values. The average difference is not 
exactly 0. This comes from the fact that the average veloc- 
ity value used for the normalization includes presampler 
data, while the gap calculation presented in Figure [Ml con- 
tains only layer 2. 



The analysis presented 
HV = 2 kV and 
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One can see that the ratio between calculated and 
design values, spanning a gap range between 1 mm and 
2.5 mm, has an RMS of 0.83%, i.e. typically 16 /im. In the 
presampler, the corresponding dispersion is 7 /im, reflect- 
ing a more rigid fixing of the electrodes defining the gaps. 
In the barrel part one recognizes the systematic effects 
in the results discussed in Section 16.1.11 (see in particular 
Figure [TU)) associated with the slight bulging of the ab- 
sorbers, and the "transition regions" at r\ = 0, ±0.8 and 
±1.4. Strictly speaking these transitions areas, for which 
additional effects enter into play, should be corrected for 
in the calculation of the RMS. In the endcaps the sta- 
tistical power is unfortunately lower giving rise to larger 
fluctuations, but no significant trend is observed. 
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Fig. 24. Relative difference between the design gap values and 
the values extracted from Tdrift measurements. 



Figure [25] shows the drift velocity obtained using Equa- 
tion [5T] as a function of pseudorapidity and the same nor- 
malization as above. As opposed to Figure [2U which gave 
the velocity at a reference field of 1 kV/mm, Figure [25] 
shows the drift velocity at the local operating field, which 
is directly related to the peak current (see Section [5]) as- 
sociated with an energy deposition. 



In the barrel region, the drift velocity is essentially flat, 
with a slight modulation reflecting the variation of the ab- 
sorber thickness with pseudorapidity. Taking the average 
value of the drift velocity in sectors of Ar\ x A<fi — 0.1x0.1, 
as was done in Section [6. 21 for Tdrift, one obtains a distri- 
bution with an RMS of 0.29% exactly equal to what was 
derived in Section 16.21 from the RMS of the T^h /t distri- 
bution, showing the expected consistency of the analyses 
using T drlft or V dr ift- 

In the endcap region, one observes the 6 sawteeth on 
each side resulting from the finite granularity of the HV 
distribution (see Figure [3]) ■ Corrections are made in the 
energy reconstruction to normalize the response of each 
strip in pseudorapidity to the response of the strip in the 
center of the HV sector, using the power law dependence. 
Beside these modulations, one observes that: 

— the average velocity in the endcaps is smaller than 
in the barrel. In the energy reconstruction this is ac- 
counted for by correction factors (which also take into 
account the fact that the lead thicknesses are differ- 
ent) determined from test beam and implemented in 
the detailed Monte Carlo simulation of the full ATLAS 
detector. 

— the measured velocity averaged over an HV sector some- 
what diminishes with increasing pseudorapidity. This 
effect goes in the same direction (lowering the response) 
as the reduced contribution of liquid argon to shower- 
ing/conversion effects at large pseudorapidities (small 
gaps). Both effects are qualitatively counterbalanced 
by the fact that the relative contribution of bends 
as compared to flat parts is lower at high pseudora- 
pidity, resulting in an increased response. As already 
mentioned, detailed Monte Carlo simulations normal- 
ized with test beam scans were used to determine the 
HV values optimizing the uniformity of response of 
the endcaps. This will be cross checked when enough 
Z° — > e + e~ decays become available. 



11 Conclusions 
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Fig. 25. Drift velocity versus r\ in the layer 2 at the operating 
point extracted from Tdrift measurements. 



We have shown in this paper that sufficient amounts of 
ionization data (~ 0.5 million pulses of energy larger than 
~ 1 GeV) can be used for a precision measurement of the 
average electron drift time in each cell of the highly gran- 
ular LAr electromagnetic calorimeter of ATLAS that has 
been readout with fast electronics, in the current mode. In 
this regime, the recorded energy is directly proportional to 
the drift velocity of ionization electrons, which is readily 
obtained from the drift time measurement. Furthermore, 
the drift velocity and thus the recorded energy are ~ 4 
times less sensitive to gap variations than the drift time. 

Taking advantage of these facts, we derived an esti- 
mate of the calorimeter non-uniformity of response due to 



gap size variations, of (0.29 



+0 o° l) % and (0.54 



+0.06>, 
-0.04; 



% re- 



spectively for the barrel and the endcaps. The other main 
contribution to the intrinsic non- uniformity of the calorim- 
eter is the dispersion of the thickness of the lead absorbers 
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which contributes 0.18% for both barrel and endcaps [H 

EI- 

The drift time is also an input needed in order to re- 
construct the signal amplitude by optimal filtering. An 
examination of the tails of the drift time distributions sin- 
gles out "transition areas" of the calorimeter, in both az- 
imuthal or pseudorapidity angle, where the electrical field 
is lower than average due to "edge effects". Some mod- 
ulations in the third layer of the barrel have also been 
observed. 

The analysis method used to derive the drift time pro- 
vides as another parameter the average absolute value of 
the amount the electrodes are off center between their two 
neighboring absorbers. The values obtained are around 
146 fim for both barrel and endcap accordion layers, and 
are substantially smaller for the presampler (66.5 /im) as 
expected from its design. 

The drift velocity, rescaled to a field of 1 kV/ mm, is 
obtained from the drift time measurements leading to an 
average of (4.61 ± 0.07) mm/ /lis. This value is compati- 
ble with previously published measurements at the same 
operating temperature of 88.5 K. 

The measurements presented in this paper illustrate 
the accuracy achieved with this method even using cos- 
mic muon data, thus demonstrating that it can be used to 
correct for the measured gap variations in order to even- 
tually reduce the constant term of the energy resolution, 
especially if the measurements are repeated with collision 
data. It is therefore important, in the quest to improve the 
energy resolution constant term, that in the future these 
measurements be done with LHC collision data. 
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